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Abstract: 

An effective Hamiltonian and equations of motion for treating both the resonant dipole- 
dipole interaction between two-level atoms and the resonant atom-field interaction are de- 
rived, which can suitably be used for studying the influence of arbitrary dispersing and 
absorbing material surroundings on these interactions. It is shown that the dipole-dipolc 
interaction and the atom-field interaction, respectively, are closely related to the real part 
and the imaginary part of the (classical) Green tensor. The theory is applied to the study 
of the transient behavior of two atoms that initially share a single excitation, with special 
emphasis on the role of the two competing processes of virtual and real photon exchange in 
the energy transfer between the atoms. To illustrate the powerfulness of the theory, specific 
results for the case of the atoms being near a dispersing and absorbing microsphere are 
presented. In particular, the regimes of weak and strong atom-field coupling are addressed. 

PACS: 42.50.Ct, 42.50.Fx, 42.60.Da, 82.20.Rp 

1 Introduction 

The progress in trapping single atoms in microresonator-type equipments has offered novel 
possibilities of studying fundamental quantum processes on the level of a very few atoms. 
Let us consider two two-level atoms and assume that one of them is initially prepared in 
the upper state. Various processes are possible, (i) The excited atom can spontaneously 
emit a real photon that leaves the resonator or is reabsorbed by the atom, (ii) The emitted 
photon is lost because of material absorption, (iii) The excited atom undergoes radiationlcss 
decay because of material absorption, (iv) The emitted photon is absorbed by the second 
atom that is initially in the lower state, i.e., the excitation energy is exchanged between the 
two atoms, (v) The excitation energy is exchanged between the atoms via virtual photon 
emission and absorption, which is commonly termed resonant dipole-dipole interaction. All 
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these processes sensitively depend on the electromagnetic-field structure which is controlled 
by the material surroundings forming the microresonator. 

In the regime of weak atom-field coupling, the mutual interaction of atoms has typically 
been described by an effective two-body potential involving atomic variables only. By an- 
alyzing a one-dimensional cavity model, it has been shown that this concept may fail to 
give a correct description of the interaction at least in the strong-coupling regime, where the 
electromagnetic field degrees of freedom can no longer be eliminated Q . In order to describe 
both the resonant dipole-dipole interaction and the resonant atom-field interaction of two 
two- level atoms in a high-Q cavity, an effective Hamiltonian has been proposed 0. Apart 
from the fact that the introduced coupling parameters are not specified and their relation to 
each other thus remains unclear (also see |3|, Q ) , material absorption cannot be taken into 
account, because of the underlying concept of mode decomposition. 



2 Hamiltonian 



In order to treat the problem more rigorously, let us start from the multipolar-coupling 
Hamiltonian for N two-level atoms [positions r^, transition frequencies oj a , transition dipolc 
moments (A— 1, 2, N)] that interact with the electromagnetic field via electric-dipole 
transitions in the presence of dispersing and absorbing bodies J5], [| : 

r r°° r°° r n 

H= d 3 r dujfiLuf1{r,uj)f(r,Lj)+^2±hujAcrAz-Y. du} d A E(r A , w) + H.c. (1) 



(2) 



(cU = d A a A + d* A a A ), where 



E(r,u;)=i A /A^ f d 3 rVei(r^)G(r,r>)f(r>). 

V 7T£0 C Z J 



Here, G(r, r',u>) is the classical Green tensor satisfying the inhomogcneous Helmholtz equa- 
tion 



• s(r, uj) — V x Vx 



G(r,r',w) = -5(r-r' 



(3) 



with e(r, w)=£r(t, u>)+iei(r, to) being the complex (Kramers-Kronig consistent) permittivity 
of the material surroundings. The fields f (r, u) and f^(r, w) are bosonic ones which play the 
role of the fundamental variables of the electromagnetic field and the medium, including a 
reservoir necessarily associated with the losses in the medium. It is not difficult to see that 
f(r,w) obeys the Heisenberg equation of motion 



f(r,w) 



;f(r,w) + — 1 



' ei(r, lS) 
hireo 



^d A G*(r A ,r,c). 



(4) 



To treat the off-resonant atom-field interaction, we proceed as follows, (i) The w-integrals 
in Eq. (Q) are decomposed into on- resonant parts (denoted by f °° dui . . .) and off-resonant 

parts (denoted by/ °°do; . . .). (ii) For f(ui) [and ff(cd)] in the off-resonant part of the third 
term on the right-hand side in Eq. (Q), the formal solution of Eq. (0) is substituted, i.e., 



f (r, u, t) = f frec (r, u, t) + ^J £ J^£. J o df d A (t>) G*(v A , r, u)e 



-iw(i-t') 



(5) 
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where ff ree (r, w, t) evolves freely, (iii) In the resulting expression, slowly varying atomic 
operators are put in front of the ^'-integrals, so that the ^'-integrals can be performed 
to (approximately) yield ^-functions which are (approximately) replaced by their principal- 
value parts (because of the off-resonance condition) . (iv) The expectation value with respect 
to the off-resonant medium-assisted field is taken, by assuming that it is in the vacuum state. 

If the differences between the atomic transition frequencies are small compared to the 
frequency scale of variation of the Green tensor, then the procedure outlined changes the 
Hamiltonian (|l|) into the effective Hamiltonian 

H e g = / d 3 r / dto hoj f^(r, q;)f (r, ui) + \huj A <T Az 
J Jo A 

i r'°° v -i 

- V U A 'A'o\a A ' -Y. duJ d AE(r A ,w) + H.c. , (6) 

A, A' A 

where the notation ^2' A A , indicates that A' . The uj a are the shifted transition frequen- 
cies, and the 5a* A' (A^ A') are the resonant dipole-dipole coupling strengths, 

&A = u A - 5 A * A , (7) 

6a*A = $2*A - 6 A'A' ( 8 ) 

Sa*a> = S A , A , +8 A *A'i ( 9 ) 
(+) = V r ^ ^ dA l mG (r A ,r A ^ )dA , 

[V - principal value]. The notation A* (A'*) means that d A (d A r) in Eq. jl0| ) has to be 
replaced with its complex conjugate (d^,). Recalling the Kramers-Kronig relation for 
the Green tensor, we may approximately rewrite Eq. (0) as 



'A* A 1 



heoc 



^d A ReG(r A ,r A >,oj A >)d A ,, (11) 



which reveals that the resonant dipole-dipole interaction is closely related to the real part 
of the Green tensor. Accordingly, from Eqs. ([lO]) and (|J) it (approximately) follows that 

5a*a = ^2 d A Re G ( r A, r A , u A )d A - 26+, A . (12) 

Note that the Green tensor can typically written as a sum of the vacuum Green tensor Gy 
and a reflection part Gr. Due to the singularity of ReGy at equal space points, Eq. ( |l2| ) 
actually applies to the reflection part only. The vacuum part can be thought of as being 
already included in lu a . 



3 Equations of motion 

Let us assume that the atoms initially share a single excitation while the medium-assisted 
electromagnetic field is in the vacuum state. In the Schrodinger picture, we then may write, 
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on omitting off-resonant terms, the state vector of the system in the form of 



A 



)4 l^)|{o}) 



du; C L (r, u, t)e-« u -°)*\L)&(T, w)|{0}) 

(13) 

(u>= 5 ^2 a &a)- Here, \U A ) is the atomic state with the Ath atom in the upper state and all 
the other atoms in the lower state, and \L) is the atomic state with all atoms in the lower 
state. Accordingly, |{0}} is the vacuum state of the rest of the system, and //(r, w)|{0}) is 
the state, where a single quantum is excited. 

Basing on the Hamiltonian (ra), it is straightforward to derive the equations of motion 
for the slowly varying probability amplitudes C A (t) and Ci,(t). By formally integrating 
the equation for Cj,(t) under the initial condition that Ci(t = 0) = 0, and substituting 
the formal solution into the equation for Ca (t) , we obtain the following system of coupled 
integrodifferential equations for the C^(t): 

C A (t)= J2 iSA-A'e^-^CA'W + Y, dt ' dojK A *A>(t,t';tu)C A '(t'), (14) 

a, JO JO 



A 
A'^A 



where 



K AA ,(t,t';w) = - 



hlT£n 



(15) 



Equations ( [l4|) and (|l5| ) reveal that the resonant atom-field interaction is closely related to 
the imaginary part of the Green tensor. Since the real part and the imaginary part of the 
Green tensor are related to each other, it is clear that the resonant dipole-dipole coupling 
strengths and the resonant atom-field coupling strengths cannot be chosen independently 
from each other in general. 



4 Two atoms near a microsphere 

To illustrate the theory, let us consider two identical atoms (A and B) near a dispersing and 
absorbing microsphere and assume that they have equivalent positions and dipole orienta- 
tions with respect to the sphere such that the relations K,a*a = K-b*b and 1Ca*b = K-b»a 
are valid, where 

tup" 

ICa*A' = fe ^' 2 d A G(r A ,r A >,u>A')dA' = -\T a *a> + i^A'A'- (16) 

Examples of Tab and Sab as functions of frequency in the vicinity of a surface-guided field 
resonance are shown in Figs, [l] and ^[ where a single-resonance Drude-Lorentz-type dielectric 
has been assumed (for details, see 0). It is seen that near a microsphere resonance the 
dipole-dipole coupling strength sensitively change with frequency. In Fig. [l], the two atoms 
are located at diametrically opposite positions (outside the sphere), which is an example of 
the interatomic distances being much larger than the wavelengths. Using Eqs. ([)]) and (|l(J) 
and approximating the microsphere resonance line by a Lorentzian, it is not difficult to prove 
that the absolute value of Sab vanishes exactly on resonance and peaks at half widths of half 
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Figure 1: Tab, Sab (dashed line), and Saa (dotted line), which almost coincides with Sab, 
in units of the free-space spontaneous decay rate To, are shown for two atoms situated near 
a dielectric microsphere of single-resonance Drude-Lorentz-type [o>t, transverse frequency; 

= 0.5 wt, plasma frequency; 7=10~ 6 wt, absorption parameter; cZ=20At, sphere diameter 
(At = 2itc/ujt)', Ata = ta — 0.5 d = Ars > 10~ 3 At, distance of the atoms from the sphere 
surface; <1a = d#, radially oriented real transition dipole moments; (a) Ata = 0.02 At; 
R = 20.04 At, interatomic distance; (b) u> = 1.05048621 wt]. For comparison, the inset in (b) 
shows the free space case. 
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Figure 2: Same as in Fig. 0(a), but for a short interatomic distance i? = 0.01 At- 



maximum. Fig. |l|(b) reveals that for interatomic distances that are large compared to the 
distances that would be required in free space strong dipole-dipole coupling can be realized, 
even for moderately small atom-sphere distances. For the example considered, the dipole- 
dipole coupling strength achieved, e.g., at an interatomic distance of about 20 wavelengths 
can be as strong as that one in free space at a distance of about 0.012 wavelengths. The 
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Figure 3: Pa($) (solid line) and P B (t) (dashed line) are shown for weak resonant atom-field 
coupling, with the data being taken from Fig. |l](a) [curves 1: u> A — 1.050485 ojt, = T BB 
= 640.848r , T AB = 640.319 r , 5 AB = 1112 r ; curves 2: Q A = 1.0504867 w T , ^aa = T BB = 
8372 r , r AB = 8371.5 r , 5 AB = 0}. 



relation between the strengths of dipole-dipole and atom-field coupling sensitively depends 
on the relative positions and dipole orientations of the atoms to each other and to the 
microsphere. By moving the atoms close to each other, the contribution of the vacuum 
part to the Green tensor increases, which allows one to vary Sab while keeping Tab almost 
unchanged (see Fig. |J). Note that for other resonance lines, Tab and Sab can switch signs 
0. Clearly, the aforegiven discussion also holds for whispering gallery resonances. 
Introducing the probability amplitudes 

C±(t) = 2-5 [C A (t) ± C B {t)\ e =F*A.*t (17) 

of the superposition states |±) = 2~ 1//2 (\Ua) ± \Ub)), from Eqs. (14) we find that the 
equations for C+(t) and C-{t) decouple, 

C±(t)= dt' duK ± (t,t';uj)e Tl5A * B( - t ' t ' ) C ± (t'), (18) 
Jo Jo 

K±(t,t';Lu) = K A « A (t,t';uj)±K A , B (t,t';uj). (19) 

The integrodifferential equation ( |T8| ) should be solved numerically in general. Closed solu- 
tions can be found in the limiting cases of weak and strong atom-field interaction. 
In the weak-coupling regime the Markov approximation applies, thus 

c±(t) = 2 - 1 / 2 e- r±t / 2 , r± = r A » A ±r A * B . (20) 

The upper-state occupation probabilities P A ( B )(t) = \C A (B){t)\ 2 then read 

P A (B)(t) = ^ [cosh (T A , B t) + (-) cos {2S A , B t)]e- rB * Bt . (21) 

From Eq. ( fjl| ) it is seen that when |<5a*b| 3> T B * B , then the excitation energy is ex- 
changed back and forth between atoms A and B (Fig. @, curves 1). In the opposite case 
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Figure 4: Pa(^) (solid line) and Pb{^) (dashed line) are shown for strong resonant atom-field 
coupling, with the data being taken from Figs. |l](a) and |^ [To = 10~ 6 ujt, f2+ ~ 128 To, Aui m 
~5x1CT 7 cjt; (i) u) A = 1.04835747 w T , i?~0.01A T , 6 AB ^-2129T ; (ii) ^ = 1.05045444^, 
i?~0.027A T , 5ab^-32.2T ; (iii) = w m = 1.0504867 w T , i?=20.04A T , ^ 0]. 



of \Sa*b\^^b-b, Pa and Pb quickly approach equal values and then decay very slowly 
(curves 2). During this slow decay process, the two atoms are partially entangled |7). Both 
types of temporal evolution can also appear in free space, but only for interatomic distances 
that are much smaller than the wavelengths. 

In the strong-coupling limit, we restrict our attention to the case when the absolute value 
of the two-atom term Ka* t'; to) is of the same order of magnitude as the absolute value 
of the single- atom term KA*A{t, t';u), so that there is a strong contrast in the magnitude of 
K + (t,t';uj) and K-(t,t'-,uj). As a consequence, the strong-coupling regime can be realized 
for either the state |+) or the state |— ), but not for both at the same time. Assuming 
that the field resonance strongly coupled to the atoms has a Lorentzian shape, with u> m 
and Auj m being the central frequency and the half width at half maximum respectively, we 
can perform the frequency integral in Eq. ( |l8| ) in a closed form, on extending it to ±oo. 
The further calculation can then be performed as described in 0. In particular for exact 
resonance, i.e., uj m =uja T^a*b, we derive 

C±(t) = 2-h- AuJmt/2 cos(fl±t/2), fl± = y/2T±Aw m . (22) 
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For the probability amplitudes of the remaining states |=p), which are weakly coupled to the 
field, we have C T (t) = 2- 1 /2 e -r T t/2 It then f u ows ^at 

PA(B){t) = i[e- r T t + e - A "'"*cos 2 (r! ± t/2) + (-) 2 e -( A ^ +I > )*/ 2 cos(fi±t/2) cos(25 A * B t) 

(23) 

The upper (lower) signs refer to the case where the state |+) (|— )) is strongly coupled to the 
medium-assisted field. Typical examples are shown in Fig. ||. When the atoms are sufficiently 
close to each other (4|5^«s| 3>^±), then a beating-type behavior ofPA(t) andPe(t) as shown 
in Figs. ||(iA) and (iB), respectively, may be observed. The beat is between the oscillation 
of the frequency 2d a* b, which arises from the two-atom dipole-dipole coupling, and the 
Rabi oscillation of frequency fl±. Recall that the collapses and revivals of the atomic level 
populations in the Jaynes-Cummings model are caused by the presence of oscillations with 
noncommensurate single-atom Rabi frequencies. Partial trapping of the excitation energy in 
atom A as shown in Fig. |^(ii) may be observed if the distance between the atoms is slightly 
increased (4|(5^4* ^ | ~f2-|-). This trapping can be understood as resulting from a destructive 
interference between the two channels of energy transfer, one via virtual and the other via 
real medium-assisted field excitation. Finally, Fig. |](iii) shows that for larger distances the 
motion becomes governed by the Rabi oscillations (4|<5yi»_B| ~0^;f2±). 

5 Outlook 

The present analysis has left a number of open questions, on which future work will con- 
centrate. Depending upon the strengths of the resonant dipole-dipole coupling and the 
resonant atom-field coupling and their relations to each other, multiplet spectra of the emit- 
ted light can be expected. Further, the problem of mutual interaction of atoms whose 
transition frequencies must be regarded as being different with regard to the variation of the 
medium-assisted Green tensor needs special emphasis, because in such a case the effective 
Hamiltonian @ does not apply. 
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